DE-FG05-92ER40717-46 



On 
On 



^1" 



Oh 

D 



X 



A<P 4 Theory from a particle-gas viewpoint 



M. Consoli 



Istituto Nazionale di Fisica Nucleare, Sezione di Catania 



q ' Corso Italia 57, 95129 Catania, Italy 



and 

P. M. Stevenson 



> 

ON . T. W. Bonner Laboratory, Physics Department 

Rice University, P.O. Box 1892, Houston, TX 77251-1892, USA 

r-> 

ON 

c | ■ Abstract: 



We discuss the physics of the 3 + 1 dimensional A$ 4 quantum field theory in terms 

ft, 

of the statistical mechanics of a gas of particles ('atoms') that interact via a — l/r 3 -plus- 
hard-core potential. The hard-core potential, 5^ 3 \r), arises from the bare vertex diagram, 
while the attractive, long-range — 1/r 3 potential is due to exchange of a particle pair via the 
t, -u-channel "fish" diagram. (Higher-order diagrams preserve this form of the interparticle 
potential.) For sufficiently small atom mass, the lowest-energy state is not the 'empty' state 
with no atoms, but a state with a non-zero density of spontaneously created atoms, Bose- 
condensed in the zero-momentum mode. This corresponds to the spontaneous-symmetry- 
breaking phase transition, and the 'phonon' excitations of the Bose condensate correspond 
to Higgs particles. The important point is that the phase transition happens while the 
atom's physical mass m is still positive: it does not wait until m? passes through zero 
and becomes negative, contrary to the assumption of a second-order transition, on which 
renormalization-group-improved perturbation theory is based. 
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1. Outline 



In this paper we consider the physics of (A<3? 4 )3+i theory from a "particle" viewpoint - in 
the same sense that quantum electrodynamics (QED) can be viewed, not as an abstract 
field theory, but as a theory of charged particles interacting via a 1/r potential. Our aim is 
to gain some really 'earthy' intuition about the physics. We begin by outlining the whole 
story, the component parts of which will be discussed in the remainder of the paper. 

It is well known that the non-relativistic (NR) limit of A$ 4 theory corresponds to 
spinless particles that interact via a 5^ 3 )(r) potential Q. Our key observation is that a 
qualitatively different potential arises from the exchange of a particle-antiparticle pair via 
the "fish" diagram. For massless exchanged particles this gives rise to an attractive, long- 
range, —1/r 3 potential. (Including the mass yields an additional factor of 2mrKi(2mr) 
which supplies an exp(— 2mr) factor at large r.) Thus, the net interparticle potential has 
a "— l/r 3 -plus-hard-core" form. The hard core, a regularized form of the <5^ 3 ^(r) term, has 
large positive potential for r < tq, where 1/ro acts as an ultraviolet cutoff. 

We therefore consider the statistical mechanics of a dilute gas of N such bosons in a 
box of volume V. Without the —1/r 3 term this would be the classic hard-sphere Bose-gas 
problem ||, |3j, |J. There is interesting physics - Bose-Einstein condensation, corresponding 
to SSB from a symmetric phase with 'atoms' to a broken phase where 'phonons' are the 
excitations. The spectrum of the phonons differs from that of the atoms. With NR 
kinematics, where the atom spectrum is k 2 /2m, the phonons have a Bogoliubov spectrum 
(k/2m)^J k 2 + 16-irna, where a is the atoms' scattering length and n is their number density, 
N/V. Modifying the analysis to allow for relativistic kinematics yields a phonon spectrum 
that, for m small, is of normal relativistic form E(k) = \Jk 2 + 8irna. The parameter 

M(n) = V8irna (1.1) 

is then the phonon mass, and we can identify the 'phonons' as Higgs particles. The non- 
trivial change of spectrum survives even in the limit a — > 0, provided that na remains 
finite H. Although this requires n — ► oo, the gas is still dilute because na 3 vanishes; i.e., 
the spacing of the atoms is vastly greater than their scattering length. In such a limit the 
approximations na 3 <C 1 and ka <C 1 required by the analysis become exact. One then has 
a "trivial" -but-not-entirely trivial theory |J|: the phonons are exactly non-interacting, but 
are non-trivial coherent states of indefinitely many atoms, and their spectrum is manifestly 
not the original atom spectrum. 

All of this assumes, however, that the system is constrained to have N particles (with 
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the interesting physics requiring n large ~ l/°0- Relativistically, though, the 'atoms' can 
mutually annihilate, so the lowest-energy state would just be an empty box. 

However, inclusion of the — 1/r 3 long-range attraction changes the situation. It pro- 
vides a potential-energy gain that could offset the rest-mass cost associated with N particles 
and the energy cost of their mutual hard-core replusion. In that case, an empty box could 
become unstable to spontaneous particle creation, leading to a new vacuum state with a 
non-zero density of 'atoms'. We find that, indeed, this phenomenon does occur once the 
atom's mass m is sufficiently small. 

The attractive potential gives a negative energy contribution |N 2 u, where u is the 
average potential energy of a pair. For a —1/r 3 potential u is logarithmically divergent, 
since one basically finds u ~ Jd 3 r(— l/r 3 )/V. Obviously, the atomic mass m would cut 
off the integral at l/(2m), but if we make m smaller and smaller then eventually a more 
important consideration is that the long-distance interaction between two particles will 
be "screened" by intervening particles. The effective r max will then involve the density 
n = N/V. Specifically, we argue that the two virtual particles being exchanged have 
to propagate through the background density of atoms and multiple collisions effectively 
convert their propagators from 'atom' to 'phonon' propagators. Thus, r max is not l/(2m) 
but l/(2M(n)), once we are in a regime where M(n) is much greater than m. 

The energy of the gas at zero temperature (when almost all particles have zero 3- 
momentum) then consists of three terms (i) rest masses, (ii) repulsion energy (obtained 
in the hard-sphere Bose-gas analysis |3], Q), and (iii) energy gain from attraction: 

E = Nm + ci— + ±N 2 u, (1.2) 
m 

where c\ is a numerical coefficient and it is 1/V times the strength of the —1/r 3 potential 
(proportional to — A 2 /m 2 ), times k^rj^/ro) ~ — ln(naro). Dividing by the volume, this 
gives an energy density for a state of density n of the form 

£ = E/y = nm + ci h c<i — ^n 2 ln(narQ), (1-3) 

m m 

where C2 is another numerical coefficient. The absence of n 3 or higher terms reflects the 
assumption that the gas is dilute (na 3 <C 1), rendering 3-body interactions negligible. 
Indeed, the last term, resulting from pairwise attractive interactions, would have been a 
simple n 2 term (absorbable into a re-definition of a) if the long-range potential had fallen 
off faster than —1/r 3 . As it is, the incipient divergence results in the Inn factor. 

Upon minimizing £ with respect to n, we see that, once m is small enough, a state with 
non-zero density, n v , will be energetically favoured over the 'empty' vacuum, n = 0. The vi- 
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tal point is that the phase transition is first-order; it happens while m 2 is still positive, and 
involves a discontinuous change in the density from to n„. The usual 'renormalization- 
group-improved perturbation theory' picture assumes a second-order phase transition oc- 
curring only when m 2 (the renormalized mass-squared of the symmetric vacuum) reaches 
zero; i.e., when the n = symmetric vacuum has tachyonic excitations, rendering it locally 
unstable. 

The above description has a direct translation into quantum-field-theory (QFT) lan- 
guage; the density n corresponds to ±m(j) 2 , where (f) is the expectation value of the field, 
and the scattering length a is related to the coupling strength by a = A/(87rm). (Hence, 
M 2 {n) = 87ran translates to \\4> 2 , as expected.) The energy density as a function of 
density translates directly into the effective potential of QFT. It consists of a ±m 2 cf) 2 term, 
a term, and a A 2 4 ln(iA</> 2 /A 2 ) term, where the ultraviolet cutoff A is proportional 
to l/r . 

The effective potential thus has the same form as the famous "one- loop" result ||. 
Conventional wisdom holds that in A<3? 4 theory the one-loop result is not reliable — and, 
indeed, from a loop-expansion perspective it cannot be justified. Nevertheless, as we 
have argued previously this form is actually effectively exact. The particle-gas 

picture clarifies the reasons for this exactness and provides new physical insight. 

In order for the Higgs mass = M{n = n v ) to be finite we need A to tend to zero 
as l/ln(A/M/ l ) ]5[,P]-]TT|]. For the phase transition to occur, the mass term nm = \m 2 (f> 2 
must not dominate the other terms; this requires m 2 < Aikf?; i.e., m is infinitesimally 
small in the physical units defined by M^. It follows that a will tend to zero, giving 
"triviality." The new ground state is still dilute (n v a is finite, but a — ► 0, so n v a 3 — > 0), 
and so the approximations involved all become exact. Moreover, since A is of order 1/lnA, 
higher-order diagrams cannot materially alter the "— l/r 3 -plus-hard-core" form of the 
interparticle potential (see Sect. 9). To incorporate all-orders effects we simply let 'A' be 
the effective strength of the short-range repulsion (such that a = X/(8irm) is the actual, 
not the Born-approximation, scattering length). 

The "particle-gas" approach gives a clue to what goes wrong with the 'renormalization- 
group-improved perturbation theory' approach. That approach assumes that loop dia- 
grams basically renormalize the strength of the <5^ 3 ^ (r) vertex, turning the bare coupling 
into a running coupling. However, that is only part of the story; the t, u-channel "fish" 
diagram produces qualitatively different physics; long-range attraction rather than short- 
range repulsion. Higher orders then serve to renormalize both those interactions. It is as 
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though there were two operators in the game, not one. 

2. Preliminaries 

We now turn to a detailed exposition of the picture just outlined. We shall consider 
the single-component A<3> 4 theory, in which there is only a discrete symmetry — ► — <£, 
and the particle is its own antiparticle. (Generalization to the 0(N) case should be 
straightforward.) The Hamiltonian is: 



where : . . . : indicates normal ordering so that the parameter m is the physical mass of the 
'atoms.' We use the term 'atoms' as shorthand for "particles of the symmetric phase" 
(i.e., elementary excitations above the ($) = vacuum). 

The aim of this paper is to gain physical insight into results previously obtained using 
field-theory formalism. What we are after is not a better calculation, but some visceral 
physical intuition. For calculational purposes the "field language" is preferable, but for 
intuitive insight the "particle language" is uniquely valuable. We basically seek the "trans- 
lation" between these two different languages. The "particle language" is obviously natural 
for describing the NR limit of the theory, and in this sense we are following the lead of 
Beg and Furlong, Huang, and Jackiw Q. However, we wish to use the particle language 
without necessarily invoking the NR limit. As always, this is somewhat problematic, since 
"relativistic quantum mechanics" is not a completely well-defined theoretical framework. 
We shall not worry overmuch about numerical coefficients, since it is the general form of 
the result that is the key point. 

To illustrate the above remarks, we briefly mention a precedent in QED. Consider the 
electron-nucleus bound-state problem: the NR limit is simple in "particle" language - 
one just needs to solve the Schrodinger equation with a — 1/r potential. To incorporate 
relativistic kinematics and spin effects one can replace the Schrodinger equation with the 
Dirac equation for a —1/r potential. However, that does not account for all relativistic ef- 
fects; it omits effects due to virtual particle-antiparticle pairs, which give rise to the Lamb 
shift. However, even those effects can be patched into the particle language; one first cal- 
culates how virtual pairs modify the — 1/r potential at short distances, and then considers 
the additional potential term as a perturbation. Such a particle-language treatment, orig- 
inally due to Bethe [L^ ], captures the essential physics, though it is hard to make it into 
a fully satisfactory calculation (see p3| , Sect. 8.7). While a field-language approach is 




(2.1) 
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preferable for achieving a systematic calculation, the particle-language approach remains 
precious for the insight it provides. 

The plan of the paper is as follows: Sect. 3 discusses the interparticle potential; Sect. 4 
reviews the hard-sphere Bose gas and its "relativization" ; and Sect. 5 considers the effect 
of including a — 1/r 3 potential. These ingredients are combined to obtain the energy- 
density expression in Sect. 6. This is not meant to be a proper calculation; rather, it 
is a collage, cut-and-pasted from well-established sources, intended to reveal the essential 
physics. We shamelessly ignore complications that only affect numerical factors, such as 
identical-particle factors of 2. 

Sect. 7 discusses the 'translations' relating the scattering length a and density n to the 
QFT quantities A and (f) 2 . The translated results then yield the QFT effective potential. 
Sect. 8 examines the phase transition, and the continuum limit (how m and A must scale 
in the cutoff — > oo limit, and what this means for n and a). Sect. 9 shows that higher-order 
diagrams do not materially alter the interparticle potential. Our conclusions and outlook 
are summarized in Sect. 10. 

Appendices A and B review quantum-mechanical (QM) scattering from a short-range 
potential, and from a — l/r 3 -plus-hard-core potential, respectively. Appendix C discusses 
field re-scaling in particle-gas language. Appendix D offers an alternative version of the 
calculation which, if less physically insightful, is tidier and has the numerical factors 
properly in place. 



3. The equivalent interparticle potential 

The equivalence between photon exchange and the 1/r Coulomb potential (and between 
pion exchange and the Yukawa potential e~ mr /r) is well known. Quite generally, for a given 
QFT matrix element M for 2 — > 2 scattering, we can ask "What interparticle potential 
V(r), when used in QM scattering theory, would yield the same scattering amplitude?" A 



sophisticated discussion of this question is given in the review of Feinberg et al [ 14 ] . For 
our purposes a rather simple-minded version will suffice. Consider the elastic scattering 
of equal mass particles, pi,P2 — ► Pi,P2> m the centre-of-mass frame (pi = — p2 = p). 
Define E = \/p 2 + m? and let q be the momentum transfer p' x — pi, and the scattering 



angle be 6. Let us compare the QFT formula for the differential cross section [15]: 

da 1 1 



dfi 256vr 2 E 2 



\M\ 2 , (3.1) 
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with the corresponding formula in QM: 

%=\m\\ < 3 - 2 > 

where f(6), the scattering amplitude, is given, in Born approximation, by 

f(0) = -Ell f d 3 re- icl - r V{r). (3.3) 



2vr 

Note that we have "relativized" the familiar NR QM result by replacing the reduced 
mass \i = m/2 by E/2. This corresponds to replacing the Schrodinger equation with 



the "relativized" version E\ + E2 + V(r), with E{ = J pf + m 2 , where pj is particle i's 
momentum operator [|l4|, |l~6| ]. 

Comparing the two formulas (and fixing the phase factor to get a positive result in 
Q3,6| ) below) gives 

Jd\e-^V{v) = ^M. (3.4) 
Taking the 3-dimensional Fourier transform of each side gives 

Note that this "equivalent interparticle potential" is energy dependent, so it is an "effec- 
tive" concept only. However, this caveat becomes unimportant in the NR regime, E ~ m. 
Following Feinberg et al fl4j we regard the 1/(AE 2 ) factor as independent of the momen- 
tum transfer q. Hence, V(r) is just a function of the relative postion r (conjugate to q) 
that also depends, parametrically, on E. 

For A<& 4 theory the bare vertex diagram gives A4 = A and the resulting interparticle 
potential is 

V(r) = ^L\5^(r). (3.6) 

We now ask what interparticle potential is produced by the t-channel "fish" diagram. 
First we note that if A4 depends only on q =| q | then the angular integrations in ( |3.5| ) 
yield 

1 1 f°° 

= ^2^2^7 dyy S myM(q = y/r). (3.8) 

Note that V(r) is spherically symmetric and naturally has a 1/r 3 factor. 

The matrix element for the i-channel "fish" diagram, neglecting the masses of the 
exchanged particles for now, is just 

M t - cxch (q) = ^Hq/A), (3.9) 
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where A is an ultraviolet cutoff. Substituting in Eq. ( |3.8| ) yields an integral that is not 
properly convergent, but which can be made well defined by including a convergence factor 
e~ ey and then taking the limit e — > 0. In this sense we have 



and 



f 

Jo 



dyy sin y = 0, 



dyy In y sin y 



IT 

2' 



(3.10) 



(3.11) 



(We made use of formulas 3.944.5 and 4.441.1 from Gradshteyn and Rhyzik (GR) |I7|] . 
Note also GR 3.761.2, 3.761.4, and 4.422.) Eq. ( 3.10| ) implies that any constant term in the 
matrix element will give no contribution to the potential at r ^ 0. Such constant terms do, 
however, lead to a <5^(r) contribution, as one sees by returning to the form (3.5). These 
delta- function contributions can be viewed as renormalizing the bare coupling strength. 
(The s-channel diagram also contributes solely to the coupling-constant renormalization.) 
Another consequence of ( 3.1C| ) is that no dependence on A survives in the — 1/r 3 potential. 
The result of substituting (|3.9|) into fljPj ) is thus: 

A 2 



V, 



t— cxch 



(r) 



1 



256vr 3 £ 2 r 3 ' 



(3.12) 



Allowing for the mass of the exchanged particles leads to [18] 



M 



t— exch 



(<?) 



A 2 



32^ 2 



/(g) In 



/ /(g) 



+ 1 



V/(<?)-i 



+ const. 



(3.13) 



where f(q) = a/1 + 4m' 2 /q' 2 and the constant term contains a In A 2 divergence, as in 
Substituting in ( |3lj| ) we obtain the result in Q3.12 ) multiplied by 



1 

7T Jo 



dy sin y yy 2 + cr 2 In 



Vy 2 + <T 2 + y 



where a = 2mr. Thanks to ( [3.10 ) this integral satisfies 

siny 



(3.14) 



da 



a 
vr Jo 



dy 



Tn 



yjy 2 + a 2 + y 



\Jy 2 + a 2 \\Jy 2 + cr 2 



(3.15) 



which can be evaluated by using GR 3.775.1 and noting that X v — Y v — * i/\n(X/Y) as 
v — * 0. The result is — aKo(a), where Kq is the modified Bessel function of order zero. 
Integrating with respect to a gives 



I = aK l (a). 



(3.16) 
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(We used formula 9.6.26 of Ref. [19]. The constant of integration must be zero to agree 



with the massless case (a = 0) evaluated earlier.) Therefore the massive case yields Eq. 



( 3.12 ) times a factor 2mrKi(2mr). This factor becomes unity when mr — > 0, while for 
large mr it tends to ^Jnmr e~ 2mr . The exponential factor makes good sense; it is just like 
the Yukawa factor except that, as there are two particles exchanged, we get e _2mr rather 
than e~ mr . 

[One can also evaluate Vt_ exc h from Feinberg et a/'s formulas: Their spectral function 
P2 is given as —I(t)(ft(s,t)/(4:Tr) (see their Eqs. (2.52)-(2.54)), where the (ft factor here 
reduces to just A 2 , and I(t) for equal masses becomes Vt — 4m 2 / (8\/i)- Inserting this into 
their Eq. (2.28b) (times m 2 /E 2 to convert their U(r) to V(r)) gives 

1 1 f°° ( A 2 \ \/t - Am 2 



Changing the integration variable to z = ryi leads to 

A 2 



= -266^V? / eW* 2 - 4m 2 r 2 e~*. (3.18) 



2mr 



The integral trivially gives unity for m = and is easily evaluated in general using GR 
3.387.6 to give 2mrK\(2mr), confirming our previous result.] 

In summary; for m small the potential consists of a 5^ (r) 'hard core' plus an attractive 
— 1/r 3 tail out to distances of order l/(2m), beyond which the potential is exponentially 
suppressed. 

4. The hard-sphere Bose gas 

In this section we temporarily ignore the —1/r 3 part of the potential and consider a gas 
of 'atoms' interacting via a short-range, repulsive potential. We follow the hard-sphere- 
Bose-gas analysis given in the textbook of K. Huang [||. (See also and the original 
references 0.) This analysis describes Bose-Einstein condensation in a dilute atomic gas; 



a phenomenon only recently observed in a tour de force of experimental technique [20|. 

Consider a box of volume V containing N identical bosons ('atoms'). A convenient 
variable is the number density, n = N/V. [Huang uses the volume per particle, v = V/N, 
which is just 1/n.] The atoms interact via an interatomic potential, assumed to be short 
range, but which may have any shape. The point is that, in quantum mechanics, low- 
energy scattering is insensitive to the shape of the potential and can be characterized by 
a single parameter, the scattering length a — which is just (minus) the low-energy limit 
of the scattering amplitude (21 
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The analysis assumes small a in two senses: 

diluteness: no? < 1, (4.1) 

'low energy': ka -C 1, (4-2) 

where k is the wavenumber (or momentum, since we set % to unity). The 'low energy' 
assumption means that the scattering is essentially pure s-wave. 

An important technical simplification, explained in detail by Huang, is that, in this 
approximation, one may replace the short-range, repulsive potential by a <5-function "pseu- 
dopotential" Q: 

V pscud (r) = ^Hr) ( 1) r, (4.3) 



m 



dr 



where the differential operator {-§^)r acts to "weed out" any 1/r singularity in the wave- 
function to which it is applied. The Hamiltonian of the system [see Huang (10.124) and 
(13.86)] is then: 

*^i:v^E^(,-^n, (4.4) 

This pseudopotential Hamiltonian is intended to be used in Born approximation only. 
This is important because, if treated as a real potential, a repulsive (5-function, of whatever 
strength, produces zero scattering in 3+1 dimensions. (See Appendix A). The normal- 
ization of the ^-function term is such that, in Born approximation, it gives a scattering 
amplitude —a, and thus mimics the full effect of the true potential. 
In the (NR) quantized-field representation the Hamiltonian is 



H = — — /d 3 rV t (r)V 2 V^(r) 
2m J 



Substituting 



+ ^ f d 3 n f d 3 r 2 ^(r 1 )V t (r 2 )^ 3 )(r 1 - r 2 )-^-[r 12 ^(ri)V(r 2 )]. (4.5) 
m J J or\2 



e ikr r + i 

V'W = H°k^-' [o k ,a^,J = £ kk ,, (4.6) 



where the operators, a^a^, create and annihilate atoms in the free-particle state of 
wavenumber k, this becomes 





H = k 2 a\ at H — aLa) n - 



2m , k mV £— ' p q 9r 
k PQ 



E„ik.r„ 
e a p+ k«q-k 

k 



(4.7) 

For the ground state of the system (and for any low-lying excited states) one expects almost 
all the atoms to be in the k = state. Hence, the occupation numbers should satisfy 
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n k=o ~ N an d Sk^o n k/N ^ !• Thus, one can effectively make the replacement ao = 
= \/N- Dropping terms relatively suppressed by 1/N leads to an effective Hamiltonian 
III: 

H eS = + — V' \(k 2 + 8vrna)at a k + 4vrna(ata 1 , + a k a_ k )j . (4.8) 

m 2m {—? L k k j 

k^o 

[The prime on the summation has to do with the (d/dr)r subtlety; see below.] 

The next step is to diagonalize H g q by a Bogoliubov transformation. That is, one 
introduces new operators fr k ,fr k that are suitable linear combinations of a k ,a k : 



a k 



(&k " a*&l k ) • (4-9) 



The & k &_ k and & k &l k terms in if e g are eliminated by choosing 



a fc = 1 + ^-2:^+2, x 2 = 2TOJo(fc) , (4.10) 

87rna 



where Cl>o(&) is the original atom spectrum, k /(2m). This yields 



where 



which gives 



tfcfr = + E ^( fc ) & k 6 k, (4-11) 

k^O 



u(k) = Mk) (jz^f) = Mk)^Jl + 2/x*, (4.12) 
w(Jfc) = cj (A:)^14- (167ma)/(2mw o (A0). (4.13) 



The new operators, fr k ,fr k , create and destroy the elementary excitations of the Bose- 
Einstein condensate — phonons — that are not single-atom excitations, but coherent 
states involving an indefinite number of atoms. The spectrum of these excitations (their 
kinetic energy as a function of their momentum) is given by to(k). 
The ground-state energy is given by 



En 



N Y' ak = N (l + 0(vW). (4.14) 



m m , , m 
k^o 



The first term arises simply from the p = q = k = term in the summation in (f4.7|). The 
second term arises in the diagonalization step when one makes the substitution fr k & k = 
&|,6 k + 1. The k-summation can be converted to an integral by X) k ~~ > V f d? k / (2ir) 3 . The 
integral would be linearly divergent but for the prime on the summation, whose effect is to 
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subtract the divergent term [23]. The result is that this second term is suppressed relative 
to the first by a factor of V na 3 . 

An important observation || is that in the limit a — > 0, with n — > oo such that na 
remains finite, the results remain non-trivial while the approximations (4.1, 4.2) become 
exact: The gas is then infinitely dilute, since na? = {na)a? — > 0, and the assumption 
fcCl/a becomes no restriction. In this limit the atoms' interaction vanishes, but non- 
trivial effects survive because there is an infinite density of them. One effect is the non-zero 
ground-state energy term in ([4.11 ). The other is that the system's excitations are phonons, 
non-trivial coherent states of indefinitely many atoms, whose spectrum ( 4.13| ) is not the 
trivial k 2 /2m atom spectrum. The phonons themselves are exactly non-interacting, so the 
resulting theory is "trivial" in the technical sense. The relevance of this "trivial" -but-not- 
entirely-trivial limit will become clear later. 

Now let us consider the effect of allowing for relativistic kinematics. For actual atoms 
(Hydrogen, Lithium, etc.) this would be pointless, because their mass m is always much 
greater than 1/a; consequently, since the analysis assumes k <C 1/a, one necessarily has 
k much, much less than m. However, we shall want to consider "atoms" whose mass is 
small compared to 1/a. Thus, even when the energy is 'low' in the sense of k <C 1/a it 
need not be small with respect to m. Thus, for us, relativistic kinematics will be relevant. 
(We shall speak of 'low energies' if k <C 1/a and 'NR energies' if k <C m.) 

A quick-and-dirty argument for how to "relativize" Huang's results is the following. 
(In Appendix D we shall describe a proper relativistic version of the calculation.) The 
original spectrum uo(k), instead of being the NR kinetic energy k 2 /2m should be the 
relativistic kinetic energy E^ — m, where E^ = Vk 2 + m 2 . (Recall that the relativistic 
and NR conventions for the zero of kinetic energy differ by m.) As seen in the previous 
section, the formula for the scattering amplitude in terms of V(r), (|3.3[ ), has a factor of 
E, while V(r) in (|3.5| ) has a 1/(AE 2 ) factor. Therefore, with respect to the NR case, we 
must scale a by a factor of (E/m)(m/E) 2 = m/E. Hence, x 2 in Eq. ( 4.1C| ) becomes 
2Ek(Ek — m)/(8ima), and Eq. (|4.13|) becomes 



^ / „ v / 8irna 
In the limit m — ► 0, where the atom spectrum is E^ = k, the phonon spectrum E\. tends 



to \Jk 2 + 87rna [24|, which has the usual form for a relativistic particle. We may then 
identify 

M 2 (n)=8irna (4.16) 
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with the mass squared of the 'phonon' excitation (which is the "Higgs particle" in QFT 
language) . 



5. Effect of the — 1/r 3 potential 

We now turn our attention to the effect of an attractive —1/r 3 potential between the 
atoms. The special property of a —1/r 3 potential is actually evident at the classical level 
and it is instructive to look at this first. Conveniently, the text by Reif [25] discusses the 
dilute classical gas using an illustrative interatomic potential of the form: 

V{r) = { °°' r<r ° (5.1) 

I -A/r s , r > r 

with an arbitrary power s. (A is a constant, which Reif writes as uo r 6-) The internal 
energy of the gas receives a contribution from pairwise interactions that is the number of 
pairs times the average energy of a pair: ^N(N — l)u ~ \H 2 u. Because the gas is dilute, 
3-body interactions, etc., are negligible, and also the motion of any pair should not be 
appreciably correlated with the motion of the other atoms, which simply provide a heat 
bath of temperature T. Then, the probability of a given separation r is proportional to 
the Boltzmann factor, yielding [25| 

_ J^re-^)V(r) 
JVJ3 re -/9V(r) ' 

where j3 = 1/(&bT), where ks is Boltzmann's constant. 

In the high-temperature (small (3) limit one might naively say that the exponentials 
go to unity, giving 

u= d \ V {r) = - — —5, (5.3) 



V J w V J ro r s - 2 

which produces a logarithmic divergence in our case, s = 3. A more satisfactory derivation 
of this result p{| is to re-express (|5.2|) as 

^lnf [d 3 re-^ r A. (5.4) 



9/3 

Adding and subtracting unity from the integrand, one has 



with 



Jd 3 re- pv(r) = V + /(/?), (5.5) 
I(P)= [d 3 r(e~ l3Vir) -l) , (5.6) 
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so that, up to C(l/V 2 ) terms, 

u =-^w (5 - 7) 

The integral J(/3) is related directly, by B 2 = — |J(/3), to the second virial coefficient, 
which represents the first deviation from the ideal-gas law in an expansion of the pressure 
in powers of the density: 

P -n + B 2 n 2 + C(n 3 ). (5.8) 



At high temperatures I{0) can be evaluated by expanding the exponential, yielding 

47T o . „ r°° dr 



47T r°° rlr 

m = --rl + ^A(3 — . (5.9) 



The first term is the excluded volume due to the hard core. Substituting in Q5.7| ) reproduces 
the naive result (|5.3|) , and we see again the divergence as s — > 3. 

Now, it might seem odd to focus on the high-temperature limit, when we are really 
interested in T = 0. However, the divergence at s = 3 will be present at any temperature; 
it is just easiest to see it in the high-temperature limit. [For any (5 the exponent f3V(r) 
becomes small as r — > oo, allowing us to expand the exponential to obtain the l/r s ~ 2 
behaviour at large r.] Also, the high-temperature limit in the classical analysis has some- 
thing in common with the T = quantum case. If we consider a potential V(r) quantum 
mechanically as a perturbation we obtain a contribution to the energy: 

IJd'njd 3 ^ ^ t (r 1 )^(r 2 )y(|r 1 - r 2 |)V(ri)^(r 2 ). (5.10) 

This is in the NR quantized field representation [Cf. Eq. (f4.5|), but with a potential V(r) 
in place of the 5-function pseudopotential] . For the zero-temperature ground state of the 
unperturbed gas we have simply tp(r) = const. = ^/N/V (normalized to Jd 5 r^^ = N). 
Thus, the above equation gives 

1 /N\ 2 r o 

' V / d 3 rV(r), (5.11) 



2 VV 

which is of the form ±N 2 u with 

u = - I if'rVir). (5. .12) 



just as in (5.3). The point is that in the unperturbed quantum ground state the particles 
are evenly distributed over the box, just as they are in the classical, high-temperature case. 
So in both cases we get an unweighted integral over V(r). This produces the logarithmic 
divergence when V(r) cx — 1/r 3 . 
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The crucial question now is: "What cuts off this divergence?" Clearly, the atom mass 
m can act as a cutoff, since the actual potential only behaves as — 1/r 3 out to distances 
of order l/(2m), due to the 2mrKi(2mr) factor. However, when m becomes very small, 
and this cutoff distance becomes enormous, another effect comes into play; namely, the 
possibility of other atoms "getting in the way" of the interaction and tending to "screen" 
it. The crucial point is that this provides a cutoff that involves the background density n. 

A first, over-simplistic approach, which nonetheless brings out the basic point that 
r max will depend on n, is to suppose that the interaction becomes "screened" as soon as 
it is more likely than not that a third particle intervenes between the two particles that 
are trying to interact. That would give an r max of order the average atom spacing, which 

-1/3 
is n * . 

A better argument is to say that in the "fish" diagram the virtual particles being 
exchanged have to propagate through the background medium. In so doing they will ex- 
perience multiple collisions with the k = background atoms. This produces a propagator 
G' that is a sum of all possible numbers of "mass insertions," where each insertion — iA 2 
corresponds to a collision with a background atom: 

G' = G + G(-iA 2 )G + G(-iA 2 )G(-iA 2 )G + ... , (5.13) 

where G = p 'il m -2 , giving G' = p -2_^' 2 +A' 2 ) • new propagator must be the phonon 

propagator, since these are the elementary excitations in the presence of the background 
density, n, and are the only things that propagate over macroscopic distances. Thus, we 
can identify m 2 + A 2 as the phonon mass-squared M 2 (n). This implies a cutoff distance 
r max ~ rather than l/(2m). This assumes that M(n) is much greater than m; 

otherwise the atom-mass cutoff l/(2m) would take precedence. Recalling that M 2 (n) ~ 
Sima for M(n) 3> m, we see that r max oc n _1 / 2 a -1 / 2 , rather than the n -1 / 3 of the previous 
paragraph. 

6. The ground-state energy density 

We are now in a position to make an estimate of the ground-state energy density of the 
gas as a function of the atom density n. From the last section, the contribution from the 
— A/r 3 interaction is |N 2 u, with 

1 /"rmax / A \ „ Aqr A 

u = ^J ro ("5) 4vrr dr = ln ( r —/ r °)> ( 6 - x ) 
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The coefficient of the — 1/r 3 interaction, A, can be inferred from Eq. ( |3.12 ). Including a 
factor of 2 to allow for both t and u amplitudes, we have A = A 2 /(1287r 3 £' 2 ). (For the 
ground state, with almost all of the particles in the k = mode, we may set E = m.) 
Setting r max equal to l/(2M(n)) = l/(2y/8irna) we have an energy contribution 

4vr A 2 




V 128vr 3 




(6.2) 



In addition, there is the ground-state energy found in the hard-sphere Bose-gas anal- 
ysis, namely N 2 -^^. (See Eq. ( 4.14| ).) This represents the energy cost of the hard-sphere 
repulsions in the quantum ground state. Finally, there is the rest-mass cost of having N 
particles (which, of course, is ignored in non-relativistic calculations). In total we have 

2irnci A 2 

E g . s , = Nm + N + Nn o = = ln(327rnar^). (6.3) 

Dividing by the volume yields the energy density: 

27rn 2 a 2 A 2 2 \ 

c K . s . = nm H h n —— — ^ ln(327rnar ). (6.4) 

m 12Qir z m z 

We see that this is a sum of n, n 2 and n 2 In n terms; these represent 1-body, 2-body, 
and "2-body-plus- medium" effects, respectively. Possible 3-body (n 3 ) or higher terms are 
negligible because of the diluteness assumption. 

Note that if the attractive potential had fallen off faster than 1/r 3 it would have 
given a simple n 2 term — which could then have been combined with the other n 2 term. 
That would correspond to replacing a, the scattering length due to the hard-core alone, 
by the scattering length produced by the whole potential. This observation suggests an 
alternative way to obtain £ g . s .: simply replace the "a" in the hard-sphere term 2 / Kn 2 a/m by 
the "effective scattering length" a e fr(r max ) of the whole potential (hard core plus —1/r 3 ), 
cutting off the infrared divergence at r max = l/2M(n). In Appendix B we show that 
a e ff can be expressed as a — m.Aln(r max /ro), with the parameter a now representing the 
scattering length due to the core and the "short-range part" of the —1/r 3 potential. (Thus, 
a can naturally be much greater than tq). The — m^41n(r max ) term in a e ff agrees precisely 
with the Born-approximation result — not surprisingly, since it arises from the "long-range 
tail" of the —1/r 3 potential, where the potential is very weak. Thus, our treatment of the 
— A/r s interaction using lowest-order perturbation theory in Sect. 5 is effectively exact, 
provided that we re-interpret a. 

Because of the n 2 In n term in <? g . s . , the system has a first-order phase transition. We 
shall analyze this in Sect. 8, but first we re-cast the result in QFT language. 
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7. Translation to QFT language 



What do the scattering length a and the density n correspond to in QFT terms? Clearly, 
the scattering length and the coupling A are related, since both reflect the strength of the 
interaction. To leading order, we may directly compare the QFT expression for the cross 
section, (|3.l|) , with the QM result da/dfl = a 2 . Including a factor of 2 for identical-particle 
reasons |p6|1 , we have 

(2A)2 - a \ (7.1) 



256tt 2 E 2 

or 



" (7-2) 



SttE 

In the ground state, where almost all atoms have k = 0, we may set E = m. The above result 
also follows by identifying the 5^ (r) potential obtained from the bare vertex diagram, Eq. 
(|3.6| ), with the (5^ 3 ^(r) pseudopotential of the hard-sphere analysis, (|4.3| ). The same a <-> A 
connection (in the context of an 0(2) theory, where the numerical factor is different) has 



been noted in Ref. [27| . 



The other relation we need is between the background density n and the strength of 
the background field </> = (<$). In a finite-volume box, where the modes are discrete, one 
may expand the field $(x, t) as: 

= E 7^ [«ke lk - X + ale-**] (7.3) 



where carries a time-dependent factor e %Ekt , with = \/k 2 + m? and 

[ak>«k] = <*k,k'- (7-4) 

The number operator is 

N = E a k fl k- (7-5) 
k 

Following Ref. |J , one argues that in the condensate state almost all the atoms are in the 
k = mode, so N = a\ao, and one may effectively regard ao (or rather a^ R = aoe* mt ) as 
the c-number vN. Then from (|7.3|) one obtains the expectation value 



and hence n = N /V is given by 

n = \mcj) 2 . (7-7) 
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With this identification, the trick of setting ao = = \/N is equivalent to shifting the 
field by a constant <\>. 

We may now consider both M 2 and £ to be functions of (j) 2 , rather than n. Using the 
translations (7.2) and (7.7) gives 

M 2 = 8irna = 8i(im^ 2 ) ( -^—) = ±A0 2 . (7.8) 

V 87rm / 



Similarly, the energy density Q6.4D becomes gives 

The energy density, expressed as a function of the field's expectation value, is of course the 
effective potential, V e s((p), in QFT language. We obviously recognise the first term as the 
usual mass term. Up to numerical factors [28], the second term is the (A/4!)^ 4 term of the 
classical potential and the last term is the zero-point energy of free-field fluctuations about 
the background field <j> (the "Trace-Log" term in the functional evaluation of the effective 
potential a la Jackiw p9| ). Recall that ro, our short-distance cutoff, is proportional to 
1/A, where A is a momentum-space ultraviolet cutoff. 



8. The phase transition region 

We now study the conditions needed for the system to be at or near the phase transition. 
(Away from this region the physics would be totally dull.) We prefer to use the QFT 
variables and then translate back to a and n. Because we have been somewhat cavalier 
about numerical coefficients let us start from: 

M 2 (<j }B ) = co|A0|, (8.1) 

and 

Ves = \m 2 <Pl + ci^% + <*2^f4 ln(c iA^|/A 2 ), (8.2) 

where cq, c\,C2 are numerical coefficients of order 1, whose specific values won't change the 
general picture. (We have added a l B' subscript to <j) to emphasize that it is the (constant 
part of the) bare field. Field re-scaling is discussed in Appendix C.) 

If V^ff has non-trivial minima at 4>b = then vb must be a solution of dV e s/d<f>B 
{). which gives the equation: 

m 2 + ci^| + c 2 -^v 2 b (ln(c iA4/A 2 ) + ~) = 0. (8.3) 
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Let vq denote the solution in the special case m = 0; it is given by 

2 _ 2A 2 _ 1/2 ( 32vr 2 ci X 



The original equation ( |8.3[) can then be re-written as: 

/(«£) = ~ C ^ V B Hv 2 b/v 2 ) = m 2 . (8.5) 

A graph of f(v B ) climbs steeply from zero, has a maximum, and thereafter decreases, 
becoming negative when v 2 B > Vq. Equating this to m 2 we see that: (i) If m 2 is positive 
and larger than the maximum value of /, there will be no solution. This corresponds to 
V e ff having a single minimum at (/> = 0. (ii) If m 2 is positive, but not too large, there 
will be two solutions for v 2 B . This corresponds to V e g having a minimum at <fi = 0, then 
a maximum (at the smaller vb root) and then a minimum (at the larger vb root), (hi) 
If m 2 is negative there is a unique solution, with vb greater than v$. Here the origin is a 
maximum of V e g and vb is a minimum. 

Case (i) is of little interest, having no SSB. Case (iii) is interesting and perfectly accept- 
able in QFT terms, but with m 2 negative (tachyonic 'atoms') our particle-gas description 
is inappropriate (or, at least, loses its intuitive appeal). We therefore focus on case (ii), 
where m 2 is positive but smaller than the maximum value of /. Since / is maximum at 
v 2 B = e~ 1 VQ, we need 

for a non-trivial minimum of V e R to exist. A stronger condition is needed for this minimum 
to be lower than the minimum at the origin; i.e., for V e s{vs) < 0: 



m 2 < c 2 ——^v B . (8.7) 



A 2 

128^ UB 



Recalling Eq. (|8.1|) and that = M(4>b = vb), we see that m 2 must be of order \M% 
or less. 

From ( |8.1| ) at 4>b = vb, we see that M% is of order Xv B , and v 2 B (for the root cor- 
responding to a minimum of V e g) lies between e~ 1 v 2 and v 2 . Thus, from ( |S.4| ) we find 
that 

32vr 2 ci\\ , n , 



Mi = A exp 



3A C2 



We want to take the cutoff A to infinity, but we want to keep M% finite, so that the 
phonons (Higgs bosons) correspond to finite-mass particles. The only way to achieve 
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these requirements is to have A tend to zero as A — > oo. In fact, just by re- arranging (|8.8| ) 
we see that A needs to behave as: 

A 2ci 1 
16^2 ~ 3^1n(A 2 /M 2 )' ( ^ 

Thus, A tends to zero as 1/ln A. The result of the previous paragraph that m? ~ AM^ 
means that m 2 is much, much smaller than M?. Since we require M? to be finite, which 
defines our physical mass scale, we need m 2 to go to zero as 1/ln A. We also see that v 2 B 
must go to infinity, like In A so that Xv B , and hence M 2 , is finite. 
This behaviour, namely 

A = 0(l/lnA), m 2 = 0(1/ In A), u| = 0(lnA), (8.10) 

(for A in units of M/J is needed if the physics is to remain interesting as A — * oo. If A were 
larger, then would go to infinity (in particular, = 0(A) for any finite A). If m 2 
were larger, then the \m 2 cf) 2 term would completely dominate V c g, and we would simply 
have a theory of non-interacting atoms. The interesting region is where one is sufficiently 
close to the phase transition, where the symmetric vacuum and the vb vacuum are in close 
competition. The unusual feature is that the elementary excitations of these vacua, atoms 
and phonons, respectively, have vastly different masses; to keep the phonon mass finite, 
we must let the atom mass be vanishingly small. 

Using the translations a = A/(87rm) and n = \m<j) 2 B , we see that 



a = 0(l/v / rnA), n v = O(vlnA), (8.11) 

where n v is the atom density in the SSB vacuum (f>g = vb- Thus, the scattering length 
a tends to zero (though it is much greater than the core size rg ~ 1/A.) This result 



reflects the "triviality" property of the A$ theory [30, 31|. However, n v a remains finite, 
because of course M? was required to be finite. Thus, the situation mentioned earlier || is 
automatically realized; the 'diluteness' and 'low energy' approximations ( |4.1| , fhl^ ) become 
exact (n v a 3 — > 0, and ka — > for any finite momentum k). The phonons are exactly 
non-interacting in this limit, so the theory is "trivial" in the technical sense. However, the 
phonons are non-trivial coherent states of the atoms, so some interesting physics is still 
present. 

Moreover, this interesting physics can show up 'experimentally' The energy density 
difference between the SSB vacuum and the = vacuum is finite in the physical units 
defined by M^. Therefore, if we heat the walls of our box to some critical temperature, 
finite in units of M^, we can expect to see a phase transition. In the high-temperature 
phase there is no Bose condensate and the excitations are the massless atoms. 
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9. All-orders considerations 



We began by considering the interparticle potential produced by the simplest diagrams: 
The tree-level vertex gave a 5^ (r) hard-core interaction, while the 1-loop "fish" diagram 
gave an additional — 1/r 3 interaction (times 2mrKi(2mr) if the mass is retained). What 
about higher-order diagrams? What contributions to the interparticle potential do they 
give? The answer is that in the region of interest, namely for a A that is 0(l/ln(A/M^,)), 
they basically give only contributions of the same type, namely <5( 3 )(r) and —1/r 3 . Thus, 
while higher orders can change the coefficients of these terms, they cannot modify the 
basic form of the interparticle potential. 

To see this, consider a general n-loop 4-point diagram, neglecting the mass m. In 
general this will be a sum of leading, sub-leading, . . . logarithms of A/q, where q is the 
momentum transfer: 



M n = A 



n+l 



A n (lnA/q) n + B n (lnA/q) n -' + ... , (9.1 



where A n ,B n , . . . are some numerical coefficients. It will be convenient here to use units 
with Mh = 1, so that InA/M/i may be abbreviated as In A, etc.. Substituting InA/g = 
In A— In q, and recalling that we are interested in the case where A is proportional to 1/ In A, 
we see that A4 consists of (i) g-independent terms with one or more 1/lnA factors; (ii) 
lng terms with two or more 1/lnA factors; (iii) In 2 q terms with three or more 1/lnA 
factors; (iv) etc. These must be inserted in the 3-dimensional Fourier-transform integral, 
( |3.5| ) or ([T^), to find their contributions to the interparticle potential, V{r). We consider 
each in turn, (i) The g-independent terms clearly produce 5^ 3 )(r) contributions, (ii) The 
lng terms will produce —1/r 3 contributions, exactly as in Sect. 3. (iii) The In 2 q terms, 
when we set q = y/r, involve In 2 y — 21ny lnr + In 2 r; when inserted in the y integration 
Q3,8| ) these give, respectively, a 1/r 3 term, a (In r)/r 3 term, and zero (for r / 0). Thus, 
the only new contribution to V(r) is a (In r)/r 3 term, which however is suppressed by one 
more power of 1/lnA than the —1/r 3 terms, (iv) Similarly, the ln p q terms in Ai n can 
generate (In 23 ™ 1 r)/r 3 terms in V(r), but these are suppressed by p — 1 powers of 1/lnA, 
relative to the —1/r 3 terms. 

Thus, in aggregate, there is no change in the "—1/r 3 + j( 3 )(r)" form of the potential, 
except at very short distances, r ~ 1/A (or extremely large distances, r ~ A/M%). The 
modification at very small distances is not unexpected; it will smooth the join between the 
hard-core and the —1/r 3 pieces of V{r). (Any modification at extremely large distances, 
r ~ A/M?, would be physically irrelevant since the interaction is effectively cut off at 
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distances of order r max ~ l/(2M h ).) 

To summarize: higher orders give contributions comparable in size to the lowest- 
order terms, but, up to negligible corrections, those contributions maintain the basic 
"— 1/r 3 + <5( 3 )(r)" form of the interparticle potential. Our previous analysis remains valid 
provided that we understand 'A' as an effective coupling strength incorporating these all- 
order effects. We need not worry about how to actually calculate this A in terms of the 
original bare coupling(s) p^j . This is exactly like the NR Bose gas, where the short- 
range potential, whatever its shape, can be parametrized by a single parameter a, the 
actual scattering length, and one does not concern oneself with calculating a from the 
original potential. Note, however, that we should not require the effective A to be a finite 
parameter; it and a have to be infinitesimally small. 

Finally, we note that higher orders should not upset the relation between the coefficient 
A of the long-range — A/r 5 potential and A. The long-range potential arises physically 
from two short-range repulsive interactions linked by quasi-free propagation of two virtual 
particles. If A represents the actual strength of the short-range interaction, then A is 
proportional to its square. 

10. Conclusions and outlook 

Since our opening section provides an outline of the whole picture, we give only the briefest 
of summaries here: The attractive —1/r 3 interparticle potential leads to an n 2 lnn term 
in the energy density, and hence to a cf) In cj) 2 term in the effective potential, and hence to 
a first-order phase transition, occurring before m? reaches zero. 

The form of effective potential obtained here agrees with our previous work |[] and 
is, of course, just what one would obtain in the one-loop approximation Q. However, 
our point is that this form is effectively exact (in the relevant case A = 0(1/ In A), as 
A — > oo). This is not because higher-loop contributions are negligible: in fact, all terms in 
the loop expansion then have the same size. However, as explained in the previous section, 
higher-loop contributions, with A = 0(1/ In A), give contributions of the same form as at 
one loop. 

Just as in the NR Bose- gas analysis, the result hinges, not on any perturbative assump- 
tion, but on the diluteness (na 3 <C 1) and 'low energy' (ka <C 1) approximations & 0], 
both of which become exact in the continuum limit. It might appear that our treatment 
of the —1/r 3 interaction's effect in Sect. 5 relied on perturbation theory: However, as 
explained in Appendix B, the Born approximation is effectively exact because the crucial 
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ln(r max ) behaviour arises from the tail of the potential, where it is very weak. An alter- 
native calculation that is directly a relativistic version of Ref. ||'s Bose-gas analysis is 
presented in Appendix D. 

The origin of the crucial </> 4 In </> 2 term can be viewed in at least two ways. In this paper 
it was obtained as the effect of a — 1 /r 3 interaction — a relativistic-quantum modification 
of the original interparticle potential. In the usual QFT calculation, however, it arises 
from the zero-point energy of the fluctuation field <£(x) — <j>. These two descriptions, in 
"particle" and "field" languages, respectively, are just two complementary ways of looking 
at the same physics. There is a parallel with the Lamb shift in QED which can be viewed 
as the effect on the energy levels of a relativistic-quantum modification of the electron- 



nucleus potential [12]. However, it can also be viewed as arising from the interaction of 
the electron with the zero-point fluctuations of the electromagnetic field [33]. 

Both in our previous QFT calculation and here we find that the theory is "trivial," in 



agreement with a large body of evidence [3C, 31|. That is, we find that scattering ampli- 
tudes (for 'atoms' or 'phonons') vanish as the cutoff goes to infinity. However, because the 
physics of the vacuum involves infinitely many particles interacting infinitesimally weakly, 
there is a finite effect on V e s- In QFT terms this effect is just the zero-point energy of 
the free fluctuation field <I>(x) — 4>. It makes perfect sense that the effective potential of 
a "trivial" theory should be just the classical potential plus the zero-point energy of the 
"trivial" fluctuations. This predicted form of the effective potential has been confirmed 



to high precision by recent lattice data |34|, [35], [36[] . 

As so often, the difficulty is not in understanding new ideas, but in freeing ourselves 
from the tyranny of old ones. The conventional view (see, eg., |], |37], |38|) has it that 
higher-loop effects change the </> 4 In <\> term in the effective potential to something like 
4> A I | In | . There is then no SSB until m? goes negative, and the transition is second 
order. This so-called "RG improvement" arises from a re-summation of the geometric 
series of "leading logs" to obtain a "renormalized coupling constant" that runs according 
to Xr{4>) = \r([J,)/[1 — 6Aj?(/i) ln(0//u)], where [i is some finite renormalization scale and 
b = 3/(167r 2 ). Absorbing the 1 into the scale of the logarithm gives \r(4>) = l/(61n(/i/(/>)), 
where /} is some finite mass scale. The effective potential in this approach is then basically 
A,r(^)c/> 4 . However, the Landau-pole problem, due to the absence of asymptotic freedom 
in perturbation theory, renders this approach unavoidably inconsistent, as we discuss in 
detail in Ref. S. Moreover, precise lattice data for the effective potential cannot be fitted 



by any version of the "RG-improved" formula [34, 35, 36]. 
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In the particle-gas picture, if one attempted to follow the "RG-improvement" program, 
one would apply leading-log re-summation to the 4-point function, thereby obtaining a 
scattering matrix element M (q) with a Landau pole. Its 3-dimensional Fourier transform 
would then not exist, and one could not obtain any sensible result for the interparticle 
potential, V(r). 

The particle-gas language provides a clue to what is missing in current applications 
of the RG method to the A$ 4 case. Those approaches assume that the main effect of 
higher-order radiative corrections is to make the coupling constant "run"; i.e., the form 
of the action remains basically the same at all scales, just with a A coupling that evolves. 
However, while some higher-order effects do indeed renormalize the strength of the bare 
interaction, there are other contributions — starting with the t, u-channel "fish" diagram 
- that produce qualitatively different physics, namely long-range attraction. This effect 
will appear in different guises in different formalisms, but its inclusion is crucial. 
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Appendix A: QM scattering from a short-range potential 



We briefly review quantum-mechanical scattering theory for a short-range repulsive po- 
tential, which is important background to our whole discussion. The scattering of two 
equal-mass particles is equivalent to the scattering of a particle of reduced mass /i = y 
from a fixed potential We consider a potential 

V(r)-{ " r<r °' (A.l) 

r > r , 




with a high 'aspect ratio' 2/jVor 2 , > 1 |4(J. Ultimately, we shall be interested in the limit 
ro — > where the interaction becomes pointlike. In this limit the potential becomes 

V(r) -> W<5 (3) (r), W = -vrr^o- (A.2) 

3 

For 'low energies' k -C l/?"o, the scattering is entirely s-wave. The cross section is then 
given by a single term of the partial-wave expansion: 



da , „,„m2 ' 



, y p_ - 1/(^)1" = ^2 e i5 °sin<5 ~, (A.3) 

where 5o is the s-wave phase shift, defined by the asymptotic behaviour of the radial 
wavef unction: 

x(r) = rt/}(r) cx sin(£;r + 5q) as r — ► cxd. (A. 4) 

A low-energy expansion of /ccot(5o is the famous "effective-range formula:" 

yfccoUo = -- + -r eS k 2 + 0{k 4 ), (A.5) 
a 2 

where a is the scattering length and r e g is the effective range. We may neglect the r e g term 
at 'low energies' fe< 1/a |40t| . At these low energies <5o is small and negative, 5o = —ka, 
giving a scattering amplitude f(9) = —a. 

For the specific potential (ATI) the Schrodinger equation for energy E = /c 2 /(2/i) yields 

, A sinh ar + B cosh ar r < rn, , 

X(r) = ^ . f n , (A.6) 
I C sin fcr + D cos fcr r > ro, 

where ^4, C, D are constants and 

a = a(k) = (2fi(V - E)) 1/2 . (A.7) 

The condition that ^>(r) is regular at r = requires x( r = 0) = 0> an d hence fixes B = 0. 
Matching conditions at r = ro give 

A sinh aro = C sin kro + D cos fero (A. 8) 



2 
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aA cosh aro = k (C cos kro — D sin kr®) . (A. 9) 

From the ratio of these two equations we obtain 

D [k cos Atq sinh org — a sin kr$ cosh qtq] 



C [k sin /cr sinh otq + a cos Atq cosh aro] 



(A.10) 



Comparing (|A.4|) and (|A.6|) we see that tan^o = D/C. Hence, from (A. 5) we have 



a = lim (- - 1 - ) = lim . (A.ll) 

fc^o V k cot <5 / k-^o \ kC J 

Hence, we obtain 

/ tanh arn \ . . _ „ s 

a = r 1 ^ , A.12 

with a = a(0) = (2/iVo) 1//2 . Since the function 1 — tanhx/x lies between and 1, we see 
that a < ro- No matter how large we make Vo, the strength of the potential, the scattering 
length can never exceed ro. Thus, one has "triviality" — zero scattering amplitude for all 
finite energies — in the limit ro — > 0. 

In particular, in the 5-function limit, V(r) — ► W8^ 3 '(r) one has aro oc \fWjr§ — > oo 
and so a = ro — > 0. The fact that a 5- function potential produces zero scattering in 3 (or 
more) dimensions reflects the "triviality" of (NR) A<& 4 theory Q. 

In Born approximation, however, one would get 

a = h Jfr*-*"^) = ^ r3 ° Vo = h w (A - 13) 

which is finite if W is finite. One may re-write the Born result as 

a = ^(ar ) 2 r (A.14) 

and recognize that perturbation theory corresponds to an expansion of the result (|A.12|) 
for small ar^. However, a 5- function potential actually entails a infinitely large aro, as 
noted above. 



Appendix B: QM scattering from a — 1/r 3 potential 

In this appendix we treat quantum- mechanical scattering from a — l/r 3 -plus-hard-core 
potential. We shall denote the scattering length for such a potential by o e fi(r max ), since 
it is necessary to cut off the potential at large distances to avoid an infrared divergence: 

{V r < r , 

-A/r 3 r < r < r max , (B.l) 

r > r max , 
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In the region ro < r < r max the Schrodinger equation for zero energy (x( r ) = rip(r)): 



can be solved exactly. By defining 



d 2 X , 2pA . v 

+ — o-X = (B.2) 



-2(^) V2 (B.3) 



and 

X(r) = w[z)/z, (B.4) 

the equation is transformed into 



whose solutions are the Bessel functions Ji(z) and Y\{z) [19]. 
The zero-energy solution for x(r) in the three regions is then 

-Asinhor r < rp, 

\(r) - ^ (FJx(z) + GTi r < r < r max , (B.6) 

+ D r > r max , 

where A, F, G, C, D are constants and a = \J2\iV§. The solution at large r is to be regarded 
as the k — > limit of C sin kr + D cos fcr, and so C is identified with /cC Recall from 
the previous appendix that the scattering length is the k —* limit of —D/(kC), which is 
—D/C. This slight cheat allows to obtain the scattering length without needing to solve 
the Schrodinger equation for a general energy. 
Matching the solutions at r max gives 

Cr max + D = —(FJ^Zm) + GY 1 (z m )) (B.7) 

Zm 

C = -^—(FJ 2 (z m ) + GY 2 (z m )), (B.8) 

^max 

where z m is z at r = r max . From the ratio of these equations we obtain a e s = —D/C as 

2 (FJiiZrJ + GYiiZm)) 



O'cS — '"max 



1 



Zm 

(FJ 2 (z m ) + GY 2 (z m ))_ 
For large r max (^> 2/j,A) we may expand the Bessel functions for z small [|l9| : 



(B.9) 



Ji(*) = \z + 0(z 2 ), = - (-- + z[ln(z/2) + 7 - I] + 0(z 3 

TT \ Z 

J 2 (z)=0(z 2 ), Y 2 (z) = ±-(-^-l + 0(z 2 In z)Y (B.10) 
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A cancellation of leading terms leaves 

a eff = r max ^ ^21n(z m /2)+2 7 + vr^ . (B.ll) 

The ratio F/G can be found by a similar matching at r = r$: 

F_ _ zotanhar Y 2 (zo) - 2ar Y 1 (zo) 
G Ioltq Ji(zq) — zq tanharoJa^o) 

For small tq {tq <C 2[iA) we may expand the Bessel functions for large z to obtain 

F _ k + tan(z - f) 

G K tan(>o-f )-l' 1 j 

where k = zo tanharo/(2aro). Note that we have not assumed that aro is small. Because 
of the tan(zo — ^f) factors, with zq large, the result is sensitive to zq modulo 2ir, and 
so, from the definition of z, depends on the precise ratio of 2\xA to r$. Thus, the ratio 
F/G can take any value, and depends sensitively on the precise details of the hard-core 
potential and how it joins on to the — A/r^ potential. 



Substituting for z m we may write the result (|B.ll) as 



a e ff = a - 2{iA In (r max /r ) , (B . 14) 

where a = 2fiA [ln(2/x^4/ro) + 27 + nF/G]. This can be compared with the result in the 
Born approximation which is 

a c s = Jd 3 rV(r) = a Bo rn - 2^41n(r max /r ), (B.15) 

where OBorn = §A i ^ / b r o ( see ( |A . 1 3|) ) . Thus, we may effectively use the Born-approximation 
result provided we replace aBorn by the free parameter a. We may think of a as the actual 
scattering length due to the core. (More precisely, it is due to the core and the "un-Born 
part" of the — 1/r 3 potential |41| .) Just as in the pseudopotential approach of ||, all the 
details of the short-range potential can be parametrized by the single parameter a. 

The fact that the Born approximation yields the correct lnr max term is natural, since 
this term arises from the large-r tail where the potential is weak. 

In the A$ 4 context the — A/r 3 potential arises from two hard-core interactions con- 
nected by particle-pair exchange. It follows that A will be proportional to a 2 . In other 
words, if the short-range interaction is parametrized by a coupling A = 87rma, there will be 
a —1/r 3 interaction of order A 2 . We can correctly treat the —1/r 3 interaction in Born ap- 
proximation even though we cannot rely on Born approximation to relate A to the original 
bare A. 
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Appendix C: Field re-scaling in 'particle-gas' language 

In the phase-transition region, our effective potential (|3,2j), written in terms of 4>b is an 
extremely flat function, since the 4>b = vacuum and the 4>b = i^s vacua are infinitely 
far apart (v B = O(lnA)), but differ in energy density only by a finite amount. It is 
therefore natural to want to re-define the scale of the horizontal axis; i.e., to define a 
"renormalized" or "re-scaled" field 4>r. 

To do this we introduced the following procedure |7|, ||, [K|: First we decompose the 
full field <&b(x) into its zero- and finite-momentum pieces: 

®b(x) =<f, B + h(x), (C.l) 

where / d 4 x h(x) = 0. Then we re-scale the zero-momentum part of the field: 

4>l = Z^l- (C2) 

The finite-momentum modes, however, are not re-scaled. The point is that for a "trivial" 
theory, as here, scattering theory and the Lehmann spectral decomposition require the 
wavefunction renormalization constant (in hs(x) = ^/Z~hhfi(x)) to be unity. However, 
there is no constraint on Z^, since there is no scattering theory for a zero-momentum 
mode — the incident particles never reach each other. For the purpose of rendering V e R 
a finite function of any Z^ of order In A would do: What determines the absolute 
normalization is the requirement that 

d(j) 2 R 

This is a standard renormalization condition, but the context here is rather unconventional. 
In this Appendix we clarify its meaning using the 'particle-gas' language. 

First, consider a slight perturbation of the symmetric vacuum state ("empty box"). 
We add a very small density n of atoms, each with zero 3-momentum. The energy density 
is now: 

£{n) = + nm + O{n 2 lnn), (C.4) 

where the first term is the energy of the unperturbed vacuum state (zero); the second 
term is the rest-mass cost of introducing N particles, divided by the volume; and the third 
term is negligible if we consider a sufficiently small density n. Thus, we obviously have 
the relation 

m. (C.5) 



= Mi. (C.3) 

^R=VR 



d£_ 

dn 



n=0 
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The equivalent in field language follows from the translation 

n = \mcj) 2 B (C6) 



and is that 



so that 



£{4>b) =VMb) = O + im 2 0l + 0{4>% In <f> 2 B ), (C.7) 



8 2 V C 



cff 



= m 2 . (C.8) 

6 S =0 

Now, let us consider a slight perturbation of the SSB vacuum (the box filled with a 
spontaneously-generated condensate). Before we perturb it, this state has a density n v 
of atoms, where n v is a (local) minimum of £ (n). From ( p.6|) we have the translation 
n v = \mv B . This vacuum state, though complicated in terms of atoms, is simple in terms 
of phonons: By definition it is just the state with no phonons. We now perturb it by 
adding a small density n' of phonons, each with negligible kinetic energy. The energy 
density of the perturbed state is then 

£{n') = £{ri = 0) + n'M h + . . . , (C.9) 

where the first term is the energy density of the unpertubed state; the second term is the 
rest-mass cost of the added phonons; and any other terms from phonon interactions are 



= M h . (CIO) 

n'=0 



negligible if n' is small enough. Thus, paralleling (C.5) we have 

d£_ 

dn' 

It is now natural to define a "phonon field" whose constant part, /, is related to the 
phonon density n' by 

n' = ±M h f, (C.ll) 

in analogy to ( |C.6| ). The "renormalized field" <f)R is simply this / plus a constant. A 

constant must be added if we want to have (f)R proportional to efts- Since, by definition, 
/ = when 4>b = vb, we need 

<t>R = f + v R (C.12) 

with 

- = tr - -4- < c - 13 > 

vb <pb zj 
Now we can eliminate / in favour of 4>r and re- write ( C.llj) as 

n' = $M h (<f> R - v R f. (C.14) 
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Hence, (|C.9|) can be re-written in field language as 

£{<t>R) = V eS (<f> R ) = Vesifa = v R ) + Wl^R - V R f + .... (C.15) 



The crucial condition, Eq. ( |C.3| ), follows directly from this. It is simply the field-language 
equivalent of "phonon-language" equation, ( |C.10| ), and just says that the phonon mass is, 
self-consistently, M^. It is also, of course, the SSB- vacuum counterpart to Eq. ( |C.8|) for 
the symmetric vacuum. 

The moral of this story is that the constant field <Pr — vr is related to phonon density 
n' in the same fashion that (ps is related to the atom density n. Note that there is a 
duality under 

atoms <-> phonons, 4>B 4>R ~ V R, n n \ Z<t> *~~* Z^, 1 ■ (C.16) 

Physically, this means that, we may use either 'atom' or 'phonon' degrees of freedom to 
describe the theory. Small excitations about the 'phonon vacuum' are easily described in 
terms of phonons, but are complicated in terms of atoms — and vice versa. Moreover, just 
as we can describe the 'phonon vacuum' as a complicated coherent state of atoms, we may 
also view the 'atom vacuum' (empty box) as a complicated coherent state of phonons. 

One can obtain various formulas for from the above considerations. There is no 
way to express just in terms of phonon properties of the phonon vacuum, nor in terms 
of atom properties of the atom vacuum. Z^ is about the connection between these two 
states, and requires some knowledge of the atom properties of the phonon vacuum, or, 
dually, the phonon properties of the atom vacuum. For instance, in field language, one 
may write 

„-i 1 d 2 V eS 



Translated to particle language this gives 

_! 2mn v d 2 £ 



z: 



Ml dn 2 



(C.17) 

B=VB 



(C.18) 



involving how the energy density varies around the phonon vacuum when we vary the 
atom density n. (There is also a 'dual' version of both these formulas.) Note that the first 
derivative jr- , vanishes, since £{n) has a minimum at n = n v . However, in terms of 

an n=n v 

phonons, when we vary n' about n' = (which corresponds to the state n = n v ), we get 
Eq. flC.lOp . It follows that we must have = 0; i.e., around the phonon vacuum, 
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a small variation in atom density causes no first-order change in phonon density (nor in 
the energy density). One can express in terms of the second derivative as: 

7 _ 1 2mn v d 2 n' 



z: 



M h dn 2 



(C.19) 



(Again, there is a 'dual' formula). 

Z^ can be written in yet more ways. Substituting (<pR — vr) 2 = {4>b — vb) 2 /Z$ into 
( C.14p , and then using <f>s = y/2n/m from ( |C.6| ) yields 

_M h {^i- ^) 2 



m n 

Evaluated at n = n v this gives 0/0, but evaluated at n = it gives 

M h n v 



m n'\ n=0 



(C.20) 



(C.21) 



which is a 'self-dual' formula involving the ratio (density of atoms in the phonon vac- 
uum) /(density of phonons in the atom vacuum). The correctness of the formula be- 
comes obvious if we use the n,(pB and translations, since it then reduces to just 
z 4> = v b/v 2 r- 



Appendix D: A tidier calculation 



Finally, to bring together all the pieces of the picture, we describe a relativistic version of 
Huang's hard-sphere Bose-gas analysis. We start with the Hamiltonian of the relativistic 
field theory, Eq. ((2.1Q. It is convenient to subtract a term /xN, where N is the number 
operator a { c a k- I n the end we shall set fx = 0. However, when comparing with the NR 
analysis we must set [i = m to convert to the NR convention that particle rest masses do 
not count in the kinetic energy — so that in the total energy of the system, H, we must 
subtract m for each particle in the system; i.e., -f^NR = H — raN. 



The calculation then proceeds as in Sect. 4. Instead of Eq. (4.6) for the NR field 
we have Eq. ( |7.3| ) for the relativistic field J42j|. Making this substitution and writing 
o-o ~ ag ~ \/N (which is equivalent to shifting the field by a constant, <f>) leads to 



with 



(m — n)n + 



An 2 
6m 2 



+ E 



°k ak ( Ek - * + A 



+ («ka-k + 4 a lk 



m 



(D.l) 
(D.2) 



31 



Some comments on this are in order: Our H was normal ordered (with respect to mass m), 
so the A-independent terms are rather obvious. The other terms come from J d 3 x : $ 4 : . 
Just as in the NR calculation, one needs to keep terms with at least two ao or aj factors; 



other terms will be relatively suppressed by 1/yN factors. That is (in the words of Ref. 
Q) we need only take into account interaction of particles in the condensate with each 
other and interaction of "excited" particles with particles in the condensate, neglecting 
interaction of "excited" particles with each other. Finally, we observe that in the NR 
case, substituting fj, = m, and E k ~ m, and E k — m ~ k 2 /(2m) this equation reproduces 
Eq. ( |4.8| ) (except for a different ground-state-energy term f43[| ) if we identify £ with 8nna. 
This implies the identification A = 8-Kma, in agreement with ( |7.2| ) above. 

We now proceed to introduce b^,b^, just as in Eq. (|4.9|) . Requiring the &k^-k terms 
to cancel gives 

a k = 1 + x 2 - xVx 2 + 2, x 2 = -E k (E k -//), (D.3) 



which has the same form as (4.10), except for the modified definition of x 2 (which has the 
form anticipated in the final paragraph of Sect. 4). We are then left (after a little algebra) 
with 



— = const. + ^ 

fc/0 



(D.4) 



and hence 

H cS - nfi = E + J2 Hk)blb k , (D.5) 



which is the analog of (4.11). The new spectrum E k is given by 



^ = ^-W 1 + ^rt (D ' 6) 



For fi = m this gives the right NR limit. With [i = one gets 



= yjEl + £ = v /k 2 + m 2 + C (D.7) 

which has the expected relativistic form. We may then identify the 'phonon' ('Higgs') 
mass squared as 

M 2 ((f)) = m 2 +£ = m 2 + \h 2 - (D.8) 

This is the same as before, except that we have been more careful and retained the m 2 term. 
However, this makes no real difference, since in the interesting region m 2 is infinitesimal 
in comparison with £ = iA</> 2 . 
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The £0 term in ( D.5 ) is 



E n = V 



[m 



An 2 
6m 2 



1 V — 



(D.9) 



where the term in square brackets (for fj, = 0) is just the classical potential, as one sees by 



substituting n 



\-rnc 



! . The last term, which arises in the step from ( p.4| ) to ( D.5| ), is the 



direct analog of the second term in ( 4.14 ) in the NR calculation. The NR version has m 
rather than E^ in the denominator, and its k summation would be linearly divergent, but 
for the id/dr) r pseudopotential subtlety. The relativistic analog, because of the extra 
1/E). factor is only logarithmically divergent. It represents a quantum contribution to the 
ground-state energy — and it is not hard to recognize it as the effect we have discussed, 
arising physically from the quantum-induced — 1/r 3 interaction between the atoms. 

Rewriting the previous equation (for fj, = 0) , the energy density of the ground state as 
a function of n, or equivalently, <fi, is: 

l u2 f d 3 k a k 
^:rr<~r h <> — — \< r I 



V eS 



E /V = i m y + 



(27r)3£V 

Substituting for a k from (D.3), and noting that xy/x 2 + 2 = AE^E^/ (X(j) 2 ) one has 



(D.10) 



.•n ■ - \m 2 4> 2 + A</) 4 + h(M) - h(m) 



1 



\<p 2 I (m), 



where 



h(M) = J 



d 3 k l~ 
(2^2 fe ' 



him) 



d 3 k 1 



(27T) 



Eu, 

3 2 



In (ml = 



d 3 k 1 

(2vr) 3 2E k ' 



(D.11) 



(D.12) 



In field language h(M) is just the zero-point energy contribution from a free field with a 
mass M((ft), while the last two terms in ( p. 11 ) are just the subtractions implied by normal 
ordering. These terms remove the quartic and quadratic divergences in Ji(M), leaving a 
logarithmic divergence. Explicitly, in the m — > limit the result is 

A 2 



1 



h{M) - h(m) - -X(t> 2 I (m) 



ln(iA(/»7A 2 



(D.13) 



4" r ~ vv " v 256vr 2 

Thus, the result for V e s is our old friend the "one-loop" result. 

The crucial point is that this form of V c s is effectively exact. Just as in the NR 
calculation the validity of the result requires only the 'diluteness' (na 3 — ► 0) and 'low- 
energy' (ka — > 0) approximations. The only subtlety is that the scattering length a should 
be the actual scattering length. The calculation of the actual a in terms of the original 
potential may be very hard, but is of no importance. In fact the actual shape of the 
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original potential is irrelevant, (provided only that it is short-range and repulsive); its 
effect is entirely parametrized by the scattering length it gives rise to. Similarly, here the 
'A' has to be understood as some effective coupling strength, incorporating the effects of 
the original QFT interaction (which could include < I >6 , <£ 8 , . . . , terms) to all orders. In 
the QFT case, if a continuum limit is to exist (i.e., if finite physics is to exist in the limit 
A — > oo), then a and A must be infinitesimally small, which corresponds to "triviality" 
of the theory. However, there is still non-trivial physics for the ground state, due to an 
infinite number of infinitesimally weakly interacting particles. 
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